Wavelet-based multifractal analysis of laser biopsy imagery 
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ABSTRACT 

In this work, we report a wavelet based multi-fractal study of images of dysplastic and neoplastic HE- stained 
human cervical tissues captured in the transmission mode when illuminated by a laser light (He-Ne 632. 8nm 
laser) . It is well known that the morphological changes occurring during the progression of diseases like cancer 
manifest in their optical properties which can be probed for differentiating the various stages of cancer. Here, 
we use the multi-resolution properties of the wavelet transform to analyze the optical changes. For this, we have 
used a novel laser imagery technique which provides us with a composite image of the absorption by the different 
cellular organelles. As the disease progresses, due to the growth of new cells, the ratio of the organelle to cellular 
volume changes manifesting in the laser imagery of such tissues. In order to develop a metric that can quantify 
the changes in such systems, we make use of the wavelet-based fluctuation analysis. The changing self- similarity 
during disease progression can be well characterized by the Hurst exponent and the scaling exponent. Due to the 
use of the Daubechies' family of wavelet kernels, we can extract polynomial trends of different orders, which help 
us characterize the underlying processes effectively. In this study, we observe that the Hurst exponent decreases 
as the cancer progresses. This measure could be relatively used to differentiate between different stages of cancer 
which could lead to the development of a novel non-invasive method for cancer detection and characterization. 
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1. INTRODUCTION 

The recent surge in the investigation of non-invasive techniques for cancer detection through fluorescence spec- 
troscopy,^^ Raman spectroscopy,'' elastic scattering spectroscopy,®'^ optical gated imaging, optical coherence 
tomography, diffuse optical tomography, polarization gated imaging, ^''^^'^ turbid medium polarimetry^^ and phase 
contrast microscopy^^ are some areas being actively pursued to understand the micro-structure variations through 
disease progression. It is well known that the elastic scattering spectrum contains rich morphological informa- 
tion about the biological tissue samples due to the inhomogeneity of the constituent organelles' sizes.®' ^^"^® 
The angular and wavelength dependence of the elastic scattering spectra have been used to analyze such subtle 
variations in the morphological changes.^' ^^'^^'^^'^^'^® 

In this work, we perform a Monochromatic Transmission Imaging (MTI) of the tissue samples which provides 
us with small angle scattering information of the biological samples. Ghosh et al.^ have analyzed the angular 
dependence of the elastic light scattering spectra in order to analyze the changing multi-fractality of the mor- 
phological structures in the tissues. Soni et al.,^^ analyzed the multi-fractality of the refractive index variation 
captured through phase contrast microscopy where the visible range of the electro-magnetic spectrum was used. 
The MTI presented in this work is a proof-of-concept study of the multi-fractality of elastic scattering in small 
forward angles through imaging where the contribution of the larger sized scatters is more pronounced. It is well 
known that the elastic light scattering spectra is the power spectrum of the refractive index variation. ^^^"^'^ Here, 
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we restrict ourselves to a single wavelength in order to observe the refractive index fluctuation at a particular 
wavelength. 

The use of the Wavelet Based Multi-Fractal De-trended Fluctuation Analysis (WBMFDFA) allows the use of 
the Multi Resolution Analysis (MRA) capability of the wavelet transforms to isolate trends of different polynomial 
orders. This is particularly helpful in the context of inhomogeneous size distribution of scatterers in a biological 
sample. This method has been used in various contexts like determining the mult i-fr act ality in light scattering 
spectra for pre-cancer detection^ and studying tissue multi- fractality through phase contrast microscopy. Here 
we apply this method to explore the possibility of differentiating between various stages of cancer. 

This article is organized as follows: In the next section (2) , we give a brief description of the sample preparation 
and the experimental setup. In the subsequent section (3), we briefly review the Fourier analysis and the 
WBMFDFA. In section (4), we present our observations and discuss the results of the analysis. We conclude 
with future directions in section (5). 

2. EXPERIMENTAL MATERIALS AND METHODS 

2.1 Sample preparation 

Hematoxylin and Eosin (HE) stained healthy and histo-pathologically graded neoplastic biopsy samples of human 
cervical tissues sliced into 4mm x 6mm (lateral), ~ 5 /im thick sections were prepared on glass slides for the 
experiment. The 12 healthy and 22 dysplastic tissue samples were obtained from G. S. V. M. Medical College 
and Hospital, Kanpur, India. HE staining involves nuclear staining by the application of hemalum (a complex 
of aluminium ions and haematoxylin) followed by the staining of eosinophilic structures by eosin Y."^^ The slide 
preparation involved standard tissue dehydration, wax embedding and sectioning under a rotary microtome. '^^ 

2.2 Experimental setup 

A schematic representation of the experimental setup is shown in Fig. (1). A 632. 8nm (output power 5mW) He- 
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Figure 1. Schematic diagram for the experimental setup. NDF represents a neutral density filter, A is an aperture, S 
is the sample, O is a 20X objective, Li and L2 are two coUimating lenses. EMCCD is an Electron-Multiplying Charge 
Coupled Device for collecting the images. 

Ne laser (Research Electro-Optics Inc., LHRR-0200,USA) masked by a reflective Neutral Density Filter (NDF) 
(Special Optics Inc., 9-1051, USA; two NDFs were used: a) with optical density 0.5, 31.62% transmittance and 
optical density 0.9 with 12.58% transmittance) was used to illuminate the sample S. An aperture was used to 
control the beam size to ~ 1mm. The transmitted light was collected through a 20X objective (Labomed LP20X 
Semi Plan Achro) and collimating lenses on an Electron-Multiplying Charge Coupled Device (Andor iXon3 897, 
with a pixel size 16^m x 16/im and image size 512 x 512). This recorded raw image was then subtracted by 
the image of a blank glass slide to remove the artifacts arising due to glass. The exposure time was kept at 
0.4s. The resulting image was then cropped for isolating the epithelium and the stroma. False color images of 
the epithelium for healthy and dysplastic samples are shown in Fig. (2). These images were then subjected to 
fluctuation analysis. 
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Figure 2. Blank subtracted images of the epithelium for (a) healthy and (b) dysplastic tissue samples. Healthy and 
dysplastic cervical tissue sample of the stroma are shown in (c) and (d) respectively. 

3. THEORY 

3.1 Image unfolding 

The image unfolding is a method to convert two dimensional data into one dimension by horizontally (vertically) 
concatenating the rows (columns) of the image. For example, a matrix of the form 



A = 



a b c 
d e f 



(1) 



can be unfolded horizontally as 

Ah^iabcdef) (2) 

and vertically as — {a d b e c fY' where is the matrix transpose. This method has been used to study the 
multi-fractal behavior of human cervical tissues earlier^^ where the multi-fractality of such tissues were analyzed 
through phase contrast microscopy. 

3.2 Spectral characterization 

The Fourier power spectrum of a signal X(^) is given by 



P(0) = X(0) 



my 



(3) 



where 4> and ^ are dual spaces of each other and for self-similar processes, the power spectrum is well known to 
follow a power law 

P(</.) ^ c/.-" (4) 

where a is called the power law exponent. This is related to the Hurst exponent hy a ~ 2H + 1, H E [0,1] 
which is a parameter used to describe the self-similarity of mono-fractal processes and is related to the fractal 
dimension through a — 2H + 1 = 3 — Df.^^ However, the ubiquitousness of multi-fractal processes in nature 
have been well studied and characterized.^'' The multi-fractality is characterized by a described by a spectrum 
of exponents instead of a single exponent. In the next subsection, we will briefly review the multi-fractal analysis 
through wavelet based fluctuation analysis. 

3.3 Wavelet Based Mult i- Fractal De-trended Fluctuation Analysis 

The question of multi-fractal signals has been studied by Stanley and his co workers extensively though the Multi- 
Fractal Dc-trcnded Fluctuation Analysis (MFDFA).''^ The Wavelet Based Multi-Fractal De-trended Fluctuation 
Analysis (WBMFDFA) proposed by Manimaran et aZ.'^^'^^ used the Multi- Resolution Analysis capability of 
the wavelet transforms to perform the de-trending of the signals. In this method, we initially make the signal 
X stationary by calculating the log-return series r{t) — \og{X{t + 1)) — \og{X{t)) and normalize it to get the 
normalized log-return series: 

r{t) {r{t)) 



(5) 



where {r(t)) is the time average of the log-return series and a is the standard deviation of r{t). Subsequently, 
we calculate the profile of the series through 



Y{t) = J2Rit'). 



(6) 



t'=0 



We use this profile series to extract the fluctuations through discrete wavelet transform. 

The fluctuation extraction involves a wavelet decomposition using the Db4 wavelet which has a support width 
of 7 and 8 filters."*^' The profile series can be decomposed as 



Y{t)= ^>.w+E E 



dijipijit), 



(7) 
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where, V'zj {t) is the mother wavelet Db4 and (pj {t) is the father wavelet such that it is orthogonal to the mother 
wavelet. The coefficients Cj{dij) are called the low pass (high pass) coefficients and capture the trend (fluctuation). 
The profile is reconstructed at a particular level j by taking only the low pass coefficients to extract the trend 
at level j. This trend is subtract from the profile to obtain the fluctuations at each scale. However, due to 
the convolution errors, these obtained fluctuations could have edge artifacts which are removed by performing 
this fluctuation extraction on the reversed proflle and taking the average.^' '^^'^"'^^ Then these fluctuations are 
subdivided in to Us non-overlapping segments such that Ug = [N/ s\ where the segment length s is related to the 
wavelet scale j by the number of fllter co-efficients for the wavelet used and is the length of the fluctuations. 
We obtain the q*^ order fluctuation function Fq{s) for g ^ as 



and for g = 



Fq{s) = exp 



1 
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(8) 



(9) 



The fluctuation function Fq{s) and the window size s are related by 



(10) 



where h{q) is called the generalized Hurst exponent."^* However, the dependence of h{q) on q does not make it 
the ideal parameter for the characterization of multi-fractality'*'* and hence a more sophisticated function called 
the singularity spectrum is required. 

The singularity spectrum /(/3) is related to the generalized Hurst exponent by the relations 

/(/?) = q{P~h{q)} + l, where, 

/3 = — t((7), and 
dq 

r{q) = qh{q) - 1. 

Here, r(g) is the multi-fractal scaling exponent and is defined by the standard partition function based formal- 
ism'^^' and and r(g) are related by a Legendre transform. A quantity 7 = max(/3) — min(/?) or the width 
of the singularity spectrum can be an important parameter for the characterization and differentiation of the 
multi-fractality of a signal. This parameter has recently been used to characterize the network properties of 
financial markets. We shall use this width of the singularity spectrum to analyze and differentiate between 
healthy and dysplastic tissue images in the following section. 



(11) 
(12) 
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4. RESULTS AND DISCUSSIONS 

We have shown the false color blank subtracted images of histo-pathologically characterized samples of healthy 
and dysplastic human cervical tissues in Fig. (2). The inhomogeneity of the tissue micro-structure due to 
the presence of various organelles with various refractive indices can be easily observed in the figures. We use 
the Fourier analysis and the WBMFDFA in order to probe the refractive index variations associated with the 
structural changes occurring in the course of disease progression. 
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Figure 3. (Color Online) Power law exponent |ci| for the epithelium for (a) healthy and (b) dysplastic tissue samples. 
The mean values for the healthy and cancerous epithelium (stroma) are |Qf| = 1.1922 ± 0.1279 (0.9721 ± 0.1692) and 
1.1695 ±0.1389 (0.9163 ±0.1745) respectively. 



The Fourier power law analysis is based on the mono-fractal hypothesis and proffers the power law exponent 
a which is related to the Hurst exponent H hy a = 2H + 1. The Fourier power in this case is a function of 
the spatial frequency. The values of the a for the healthy and dysplastic epithelium and stroma are shown in 
Fig. (3(a)) and (3(b)). We observe that the mean absolute a for the epithelium (stroma) is 1.1922 ± 0.1279 
(0.9721 ± 0.1692) for the healthy case while it is 1.1695 ± 0.1389 (0.9163 ± 0.1745) for the dysplatic case. This 
power law behavior in the spatial frequency domain can be attributed to the inhomogeneous size distribution of 
the scatterers in the tissue micro-structure. However, this power law exponent is not independent of the scale 
as has been observed in the case of the phase contrast microscopy^^ and implies the multi-fractal nature of the 
scatterer composition of the tissues. 



In terms of distinguishing different stages of cancer, as compared to the resuhs obtained in the analysis of the 
phase contrast microscopic images (where the whole visible spectral region (400nm-800nm) is used as opposed to 
this method where a monochromatic image at 632. 8nm is used), the mean value* of a is lower. Nevertheless, the 
trend of a for dysplastic samples being lower than that of the healthy samples is corroborated at the 632. 8nm 
wavelength. As expected, the a for the epithelium is higher than that of the stroma as- The densely packed 
structure of the connective fibers in the stroma as compared to the epithelium causes ols < a^. K comparison of 
the mean values of as and a^ for health and dysplastic tissues is given in table (1). To further verify this and to 
analyze the multi-fractal nature of the tissues, we proceed to study the images through the WBMFDFA. 
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Figure 4. (Color Online) Hurst exponents for (a) healthy and (b) dysplastic tissue samples. The mean values for healthy 
and cancerous epithelium (stroma) are K = 0.3279 ± 0.0275 (0.2152 ± 0.0384) and 0.3154 ± 0.0394 (0.2003 ± 0.0260) 
respectively. 



Figure 4 shows the calculated values of the Hurst exponent for healthy and dysplastic tissues in Fig. (4(a)) 
and (4(b)) respectively. The Hurst exponent for the epithelium Hf. and that for the stroma Hg are also compared. 
We can see that while the healthy tissues show a mean Hf. of 0.3279 ± 0.0275, the dysplastic tissues display a 
mean value of 0.3154 ± 0.0394. Similarly, the values of mean Hs are 0.2152 ± 0.0384 and 0.2003 ± 0.0260 for 
healthy and dysplastic stroma respectively. This is in agreement with the results of the Fourier analysis based 
power law calculations. As observed earlier, the iJg > Hs- The comparative values of Hs and are given in 
table (1). 

We mentioned in the earlier section that multi-fractal signals require the singularity spectrum for character- 
ization and that the width of the singularity spectrum 7 is a parameter used to describe multi-fractality of the 
data under consideration. In Fig. (5), the 7s and 7e for healthy 5(a) and dysplastic 5(b) tissues is shown. The 
mean 7e and 7s are found to be 1.0493 ± 0.7226 and 0.6257 ± 0.5827 respectively. As compared to the differences 
between a^, and as and H^, and lls\ the difference between 7e and 7s is more pronounced. Still, the trend of 
7e > 7s is followed, also implying the higher multi-fractality of the epithelium as compared to the stroma. 



Parameter 


Epithelium 


Stroma 


Healthy Dysplastic 


Healthy 


Dysplastic 


[a] 


1.1922 ± 0.1279 1.1695 ±0.1389 


0.9721 ±0.1692 


0.9163 ±0.1745 




0.3279 ± 0.0275 0.3154 ± 0.0394 


0.2152 ±0.0384 


0.2003 ±0.0260 


[7] 


1.0493 ± 0.7226 0.6257 ± 0.5827 


0.4486 ± 0.3440 


0.4095 ±0.3416 



Table 1. Comparison of the mean parameters calculated through the Fourier and the WBMFDFA. [• • •] represents the 
ensemble average over the different samples of healthy and dysplastic tissues. The parameters have been averaged over 
12 and 22 samples of healthy and dysplastic cervical tissues respectively. 



*The mean value here denotes the ensemble average or the average over all the samples. 
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Figure 5. (Color Online) Width of the singularity spectrum 7 for (a) healthy and (b) dysplastic tissue samples. The 
mean values for healthy and cancerous epithelium (stroma) are 7 = 1.0493 ± 0.7226 (0.4486 ± 0.3440) and 0.6257 ± 0.5827 
(0.4095 ± 0.3416) respectively. 



5. CONCLUSION 

In conclusion, we have compared the normal and dysplastic human cervical epithelium and stroma in an effort to 
probe the differences between their multi-fractality towards developing a possible non-invasive optical technique 
for cancer detection. The use of monochromatic imagery for this purpose provides us information about the 
micro-structural changes in the tissues associated with disease progression. 

We observe that the multi-fractality of tissues decreases with progression of cancer. Though the a or 7 are not 
drastically different, the trend of decreasing multi-fractality through the progression is sufficient to quantify the 
stage of cancer. For example, the decreasing a shows that the as the cancer progresses, the tissue micro-structure 
goes from 1// towards a flatter power spectrum. Similarly, the Hurst exponents show that with the increase in 
the cancer grades, the structural organization goes from a seemingly more random behavior to a spatially long 
term correlated behavior. The decreasing 7 implies the decrease in multi-fractality with the disease progression. 

This work is a proof of concept about the exploration of monochromatic laser imagery to understand the 
variations of the tissue micro-structure when observed at a particular wavelength as compared to phase contrast 
microscopy where the whole visible range is used. However, we believe that due to the presence of various 
fluorescing enzymes in the cancerous tissues like NAD"*", the response of the tissues to different wavelengths 
would be different and a more extensive study of the same is being undertaken which will be reported soon. A 
notable feature of this study is its candidature for in-vivo examination and characterization of the MTI which is 
also under investigation currently. 
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